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Abstract 

We investigate the consistency requirements of the next-to leading BFKL equa¬ 
tion with the renormalization group, with particular emphasis on running coupling 
effects and NL anomalous dimensions. We show that, despite some model depen¬ 
dence of the bare hard Pomeron, such consistency holds at leading twist level, 
provided the effective variable a<j(t) log(l/x) is not too large. We give a unified 
view of resummation formulas for coefficient functions and anomalous dimensions 
in the Qo-scheme and we discuss in detail the new one for the qq contributions to 
the gluon channel. 


PACS 12.38.Cy 


1 Introduction 


The understanding of small-x scaling violations in QCD has become increasingly impor¬ 
tant with the discovery of the steep rise of structure functions at HERA[|l|. Such rise 
may be related to the hard Pomeron]^ small-a: behaviour p], but may be also under¬ 
stood through large scaling violations 0 and may also £t in a hxed order perturbative 
approach p. 

The above ambiguities are partly unavoidable, because of the lack of sufficiently precise 
1^, independent measurements of small-x parton densities]^], but are partly due to the 
lack of complete understanding of scaling violations, at theoretical level. In fact, in the 
small-x region, the QCD perturbative series is affected by large logarithms - similarly to 
what happens in semihard processes [|^ - and becomes dependent on the effective coupling 
constant as{t)\ogl/x or as(t)/uj, where t = log(Q^/A^) and a; = A^ — 1 is the moment 
index. 


Therefore, a hierarchy of perturbative terms is defined in the usual way 

'y{as,uj) = 7 l + NNL - terms (1.1) 

for both the singlet anomalous dimension matrix and the coefficient functions. Further¬ 
more, since gluons couple to electro-weak probes only through quarks, we actually need 
leading (L) coefficient functions and up to next-to-leading (NL) anomalous dimensions, 
in order to have a factorization scheme independent result for DIS. 


This theoretical program has been partly achieved at quark level [8-10], but is not yet 
complete at gluon level, despite various efforts [11-14]. Nevertheless, such completion is 
needed in order to well understand scaling violations, and then be able to look at the 
small-x behaviour of initial data[^ so as to hnd the hard Pomeron rise, if any. 


In this paper we wish to present a resummation method which is valid for all the 
relevant quantities up to NL level, and to apply it to coefficients and anomalous dimen¬ 
sions in the so called Qo"Scheme []^, including the new qq contributions to the gluon 
anomalous dimension, already reported in a short note [M. The method is based on the 
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idea of k-factorization P, which is just a way of exploiting the consistency of high- 
energy factorized quantities, which are k-dependent, with Renormalization Group (R.G.) 
factorization, which is not. The all-order resummation is obtained by k-integrations of 
low order k-dependent kernels. 

Therefore, we shall proceed in two steps. Our hrst goal is to explicitate the NL 
consistency requirements between high-energy and collinear factorization with running 
coupling referring, for dehniteness, to hard processes of DIS type, characterized by a 
probe with a large scale Q and by a target with scale Qq -C Q. This will allow us to 
obtain well dehned resummation formulas in terms of kernels of BFKL type p| , up to NL 
level. 

As a second step, we shall proceed to the study of the NL kernel itself, of its eigenvalues, 
and of the resummed formulas for the qq contributions. This study is based on the phase- 
space integration of the squared QQ production matrix elements [^, combined with the 
virtual corrections available in the literature [^, Q. 

To start illustrating the hrst step, let us notice that k-factorization is due to the 
exchange of an off-shell gluon of virtuality t ~ —k^, which behaves, at high energies, 
as a Regge pole yielding a quasi-const ant cross-section. On the other hand, collinear 
factorization is due to nearly on-shell quarks and gluons, providing logarithmic scaling 
violations. 


Therefore, the two kinds of factorization are not related a priori. In unbroken gauge 
theories like QGD, a relationship can be established because the (Regge) gluon polariza¬ 
tion coincides, on shell, with one of the collinear polarizations [|^]. For instance, at one 
loop level, Regge behaviour is established, in the gluon channel, by the singular behaviour 


of the DGLAP splitting functions 


Pggi^) - 


2Ca 


Pgg{z) - 


2C, 


( 1 . 2 ) 


or in moment index space a; = iV — 1, by a correspondingly singular anomalous dimension 
matrix 
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+ 7ab + 


(1.3) 


(J 

_ CYsCa 

a,b = q,g, as = - 

TT 

However, establishing such relationship at next to leading level requires the treatment 
of several subtle points. Firstly, k-factorization is useful in the gluon channel only, while 
collinear factorization involves both gluons and quark-sea. How to recover from the former 
the two-channel picture of the latter? 

Furthermore, even in the gluon channel, the high energy exchange does not yield ah 
collinear singular gluon contributions, some of which occur at low energies, as is apparent 
from Eq. (1.3). How to relate the two gluon dehnitions and to get all subleading terms? 

Finally, the renormalization group evolution with running coupling has been always 
assumed in order to extract physical predictions from k-factorization. Is this behaviour 
really consistent with the BFKL equation and its ’’hard Pomeron” singularity at NL level? 

In Sec. 2 we shah show, on the basis of the work of Fadin, Lipatov and collabora¬ 
tors [11-14] and of our group [8-10,16], that a detailed comparison of high-energy and 
renormalization group is nevertheless possible and provides resummation formulas for the 
NL anomalous dimension matrix (and for the leading coefficient functions) in the Qo 
factorization scheme. 




The basic point is that, as it appears from Refs. [8-16] and Ref. [Iffi], s-channel 
iteration is not yet relevant at NL logs level, while t-channel iteration is provided by 
a generalized high energy integral equation, in which low energy k-dependent kernels 
occur. This result, with some modihcations, forms the basis for the representation of the 
two-scale problems with anomalous dimensions that we provide here. 


Some of the questions asked above do have a tricky answer, however. Firstly, even at 

leading twist level, we hnd that the renormalization group representation with running 

02 

coupling is valid, for a; -C 1, only in the regime of t = log ^ S> fo > with 


bujt > c 


(1.4) 
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where the critical constant c is related to the BFKL kernel. This means that we have 
to stay away from a t-dependent singularity in the uj plane, or in other words, that the 
effective variable as(t) log(l/a:) should not be too large. 

Furthermore, the bare hard Pomeron singularity, which occurs in the coefficient func¬ 
tion, and ultimately dominates the small-a: behaviour, turns out to be dependent on the 
large distance behaviour of the running coupling, due to a diffusion process towards low 
values of k^. Fortunately, this fact does not affects the R. G. factorization, valid in the 
regime (1.4) PU| . 

Finally, within the anomalous dimension representation, the two-channel picture ob¬ 
tains because of the presence of two distinct anomalous dimension eigenvalues 7 +, 7 _, 7 + 
being the leading one. While the NL hierarchy is straightforward for the 7 + contribu¬ 
tions, it requires careful handling of a whole chain of collinear singular kernels for the 7 _ 
contribution, as we shall see. 


In other words, comparing the high-energy and R. G. approaches is not straightforward 
at NL level, mostly because they resum the perturbative series in different orders. 


The contents of the paper are as follows. In the next section we set up the NL 
factorization framework with scales Q and Qq and we discuss the structure of its solution 
in the anomalous dimension regime dehned before, by singling out running coupling effects, 
and by discussing the role of the hard Pomeron singularity and its model dependence 


In the following sections we concentrate on the calculation of the anomalous dimensions 
themselves. In Sec. 3 we provide the qq emission kernel for massive quarks and in Sec. 4 we 
combine it with virtual corrections in the massless limit to get the (regularized) complete 
eigenvalue and the ensuing resummation formulas. We discuss our results in Sec. 5 and 
we provide the details of the phase-space integration and of the kernel diagonalization in 
Appendices A-G. 
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2 Next-to-Leading Equation and Anomalous Dimen¬ 
sion Regime 

2.1 Gluon density and bare Pomeron with running coupling 

Let us formally define the unintegrated gluon density by the high-energy limit of 

the six-point function in Fig. 1. Such density is a function of the moment index u = N — 1 
and of the momentum transfer of the exchanged gluon, and will eventually be measured 
by coupling it to physical probes, as done later on for the DIS structure functions at scale 

Q. 

We shall assume that the gluon density satishes, at NL level, the BFKL-type equation 
of Fig. 1, whose kernel will be written in the form 

A-„(k, k') = (k, k') = (A-„(k, k') + a.A-i(k. k')), (2.1) 

OJ OJ 

a,ik^)CA _ 1 _ 1 

TT fecu log(k2/A2) hut 

where we have factored out the running coupling at the ’’upper” scale k^, i.e., 

tt,(k2) = a,(/i^)-6a^(/i^)log^ + .... =|^61og^^ . (2.2) 

Furthermore Kq represents the leading kernel 

kO = - = - - <5^(q) / ^, (q = k - k'), (2.3) 

q j? q JA2 q 

and Ki is the NL one, which can be assumed to be scale invariant and contains no high- 
energy gluon exchange contributions. Its specihc form has been discussed at length in the 
literature [11-14,16,19], and its iV/-dependent part will be explicitated in Sec. 3. 

Because of scale invariance, we can dehne the kernel eigenvalue x{l) by the equation 
j K{k,k'){kr-^-- = x(7)(k2)"-b (2.4) 

The explicit form of y is known at leading level, yielding 

Xo(7) = 2^/’(l)-^/’(y)-^(1- 7 ), (2.5) 
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and the NL qq contribution will be given in Sec. 4, where we also prove the 

factorization of the logarithm in Eq. (2.2). 

Let us stress the point that factorizing the running coupling at scale is not neces¬ 
sarily an indication of the ’’natural” scale of the problem, which rather seems to be q^, 
the emitted gluon transverse momentum (see Sec. 3). In fact, at NL level, a change of 
scale in the factorized logarithm of Eq. (2.2) can be compensated by a corresponding 
change of the scale-invariant kernel Ki. The choice of allows us to deal with a type of 
equation already studied in the literature. 


Indeed, it was pointed out long ago |^, ^ that by writing the master equation 

71 


(k) = (k) + f A'(k,k')—J'»^(k') 

UJ I 71 


( 2 . 6 ) 


in the y-moments representation 


we get a differential equation in 7 , with the formal solution 
/ 7 ( 7 ) =exp (-Lvb-)) Hd',’ (-Afcp(y)) exp 


(2.7) 




7-^(7') 

buj 


( 2 , 8 ) 


Here we have denoted by X(y) the primitive of the eigenvalue function, i.e., 

V( 7 )=/ dix(i),Xa('l) ^ f o ! 7 ' xo ( 7 ') = 2 i/>(l) ( 7 -x) - log- 


r(l-7) 


(2.9) 

and we have denoted by fcf,"^(y) the y-moment of the (low-energy) inhomogeneous term 
in Eq. (2.6). 

Actually, the formal solution (2.8) has to be modihed in order to deal with the running 
coupling (2.2) in the small k^ region around the Landau pole. Cutting off [22-24] or 
smoothing out such a pole results in modifying the large y' region in the representation 
(2.8) by adding to the inhomogeneous term some a;-dependent contribution which contains 
the Pomeron singularity ig. 
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In Ref. we have investigated the Green’s function of Eq. (2.6) in the regime 
t S> to = log(Qo/^^)! where Qq represents the scale of the inhomogeneous term and we 
have shown that it factorizes in the form 


•^< 2 ^ = i^)9^iQo) + higher twists, (t > to) 


where 





and 







Pomeron terms. 


d'y' 



+ 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


We have also found that the position and the structure of the Pomeron singularity 
in Eq. (2.12) is dependent on the smoothing out procedure of as{t), and decouples only 
in the limit of a hard initial scale (to S> 1 ), being suppressed, in that case, by (roughly) 
powers of Is?'/Q q. 

The results above show that, in the limit of t S> to 3> 1, we have exact consistency 
with the nazve perturbative result in Eq. (2.8), because we can set 7 ~ 0 in the second 
integral of this equation. If instead to is of order unity, the structure of the large t solution 
remains factorized, but the coefficient in front is not really calculable, being dependent 
on the low t behaviour of h e- on soft hadronic interactions. 


In either case the NL BFKL equation remains consistent with the R. G.. In fact, in 
the anomalous dimension regime, defined by 

bujt > X , (2.13) 

the large t behaviour of Eq. (2.11) is dominated by the perturbative branch 'j{as{t)) of 
the saddle point equation 


1 = 



( 2 . 14 ) 
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By integrating in 7 over the spread |( 57 | ~ y—huj/x'{l) on the imaginary axis, we hnd 
the expression 

^ expF^ [7(a,(t),t)], (2.15) 

V-X'(7) 

where the exponent 

F^{x,t) =xt-^X{x) = j x{(^s{r))dT (2.16) 

is the one predicted by the R. G.. It is clear that if the limitation (2.14) is not satished, 
the representation (2.15) breaks down, due to large 7 -fluctuations. 

Then, by introducing the integrated gluon density 

= / —^^(k) = / ^e^V.^( 7 ), (2.17) 

Jo tt 2mx 

we dually obtain , by Eq. (2.16), the R. G. representation 

Gi^{Q) = H ( 7 (as(t))) exp [F^ ( 7 (a,(t), t))] K^{uj), (2.18) 


where the 1 . h. coefficient is given by 

1 1 

H{x) =- 


(2.19) 


while the r. h. coefficient K^{uj) contains the Pomeron singularity and is thus, strictly 
speaking, not calculable. 


Several remarks are in order. First, in the case several saddle points are present, a 
sum over them is understood. At NL level, we need at least two of them corresponding 
to the two eigenvalues 7 +, 7 _ (Sec. 2 . 2 ). 

Secondly, in the limit Ql 3> A^, the nafve expression (2.8) applies, and we can set 

'^''‘'^^ = ^1 exp (Txm) , (2.20) 


SO that a perturbative evaluation is possible, as in the following examples. 



Initial gluon with fixed virtuality. In this case, we use k-factorization at the lower 
scale also, so that we can hx the gluon virtuality in a gauge invariant way, by setting 
= Qo > "^ke inhomogeneous term becomes a ^-function in the k-space, or 

k{fi) = ( 2 . 21 ) 

in the 7 -space. For sufficiently large to (hujto 3> x(l/2)) we can use a saddle point located 
at 7 (as(to)) 3> ^{as{t)) at the lower scale also with 7 -spread |( 57 | = \J■, so as to get 
the approximate expression 

K%Qo) ^ -^=^=== exp [-F^ (7(a,(to), to)] (2.22) 

V-X (7(«4^o)) 

which is t-independent as expected. 

Initial probe at scale Qq. This scale can be due, for instance, to heavy quark production 
{Qo = M), or to a virtual electroweak probe. In either case one should use k-factorization 
to dehne as an off-shell cross section along the lines of Refs. [8-10]. In all these cases, 
the function has, by scale invariance, a to-dependence of the type ( 2 . 21 ) apart from 

a 7 -dependent factor, computed in some cases ||, |^, | 2 ^, which is analytic in the strip 
|/?e 7 I < 1. The argument about f-scale dependence of K goes through as before. 

Note, that the perturbative evaluation (2.22) is valid only if the Pomeron term in k^ 
can be neglected. This implies that as{to) log(l/a;) should not be too large at the lower 
scale also. The latter assumption is not required in the more general expression (2.18). 

2.2 Quark sea and anomalous dimension matrix in Qo"Scheme 

It is convenient to introduce the quark-sea contribution by the flavour-singlet part of the 
F 2 structure function itself (DIS scheme for the quark). This one can be divided in a 
direct low-energy contribution Go, which contains collinear singularities but not a; = 0 
singularities, and in a gluon exchange part (Fig. 2). The latter is evaluated by using 
k-factorization, in terms of the unintegrated gluon density dehned before, as follows 

G7(0) - Gt(Q) £ki.AQ,k,m)T‘-^(h), (2.23) 
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where 02 is the gauge-invariant off-shell cross-section introduced in Ref. [§] and m is the 
quark mass, that will eventually be set to zero. 


By introducing the 7 -representation of Eq. (2.7), the r.h.s. of Eq. (2.23) takes the 
form 

1-7-71 


N 


f 




d'ji 

2ni 






^.( 7 , 71 ), 


(2.24) 


where we have introduced the double 7 -moments of (J 2 , that is , the abelian contribution 
to the qq i7-function, already computed in Ref. for heavy flavour production. In the 
a; = 0 limit, we have 


72 ) =—r(i - 71 - 72)r(7i)r(72) 


a 


TT 


B(l - 7 i, 1 - 7 i)R(1 - 72 ,1 - 72 ) 


B{1 - 7i, 2 - 7 i)R(1 - 72, 2 - 72 ) 

(3-270(3-272) 


(1-7 ( 1 - 70 ( 1 - 72 )) 


(2.25) 


The massless quark limit on Eq. (2.24) singles out the pole at 71 = 1 — 7 of the 
expression (2.25), i.e. 

7 / 2 ( 7 )= lim i7“07i,7)(l - 7i - 7 ) = 

71^1-7 

as f TT Y 2 + 37 ( 1 - 7 ) 

=- - cos 7 r 7 7 -r--TT-r 

2n \sin 7 r 7 / (3 — 27 ) (1 — 27 ) (1 + 27 ) 

so that we hnally get, for massless quarks, in the small u limit 
G'ifi.Q) -Nf j |^i 72 ( 7 )/^( 7 )e^* + low energy terms 
which becomes, in the anomalous dimension regime, 

G7(Q) = 5^7ff2(7)G“7Q). (2.28) 


(2.26) 


(2.27) 


The saddle points 'j{as{t)) are given here by all solutions of the equation (2.14), i.e., 
as{t) 


1 = 


OJ 


[Xo( 7(«7)) + «.Xi(7(«.))] 


(2.29) 


where we have explicitly shown the NL term of the BFKL eigenvalue. 
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While at leading level we only have the contribution of the BFKL anomalous dimension 


IL 


(2.30) 


1 = 

IL = f + 2C(3) (f)“ + ... 

at NL level we expect saddle points located at the two eigenvalues of the anomalous 
dimension matrix 7 ^^, i-e., 


Cp Cp ( 1 ) 

7+ ^ Igg + -^Iqg. 7- ^ Igg " -^Igg = —^%g + ' ' ' (2-31) 

where we have used the fact that, at leading level 7^5 ~ ^'Jgg, and that the NL 'jqg starts 
at one loop level ( 7 gP). 


Here a subtle point arises. The NL hierarchy of the BFKL kernel really works only 
around the leading eigenvalue 7 +, but not around 7 _. In fact, subleading BFKL kernels of 
order n > 1 (the leading one corresponding to n = 0 ) show multiple collinear singularities 
of order ^ . Evaluating the latter around ■jp — yields genuinely subleading 

terms of relative order u"'. Instead, for small values of 7 ~ 7 _ ~ all subleading kernels 
become of the same order! 


This means that, in order to recover the correct collinear behaviour of the partonic 
Green’s functions around 7 = 7 _, we need to resum a whole series of collinear singular 
BFKL kernels. For instance, at one loop level, the analogue of /^“( 7 ) in Eq. (2.7) 
(a, b = q,g) is simply found from the DGLAP equations. Keeping frozen for simplicity 
we obtain 


fJil) = 


7 


1 - 


7 


-1 


1 - 


7+ Igg 
7 


''(i-7)(i-7) V 7 


which satishes a BFKL-type equation of the form (2.6), with 
-1 


(2.32) 


k = I 7 — F 


K = kPp, 


(2.33) 


where Pp and F are defined by the expression 


a, 

r = —Pp + F, 

OJ 


(2.34) 
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which is the NL decomposition of the one-loop flavour singlet anomalous dimension matrix 
(Cfr. Eq. (1.3)). It is then easy to check that the collinear resummation in the kernels of 
Eq. (2.32) is essential to provide the correct singular behaviour at 7 = 7 _, thus cancelling 
a variety of singularities around 7 = 0 that finite order expansions have. 


On the other hand, since we are working at NL level, the one-loop result is all we need 
to describe 7 _, which does not have enhanced higher order contributions. This means 
that Eqs. (2.32)-(2.34) are enough to find the correct behaviour around 7 _ if running 
coupling effects are introduced according to the renormalization group. 


The output of this discussion is that, in Eqs. (2.18) and (2.26) we evaluate the saddle 
point 7 + by the k-factorization formulae given before, and we evaluate 7 _ by the known 
one-loop collinear behaviour, because higher order enhancement effects are not present. 
The resulting quark-sea and gluon densities are then given by 




Gf + 



(2.35) 


where 


Gf = 

G9A ^ 


7+\/“N(7+(«db)) 


exp 


Io'y+Mt))dt 


exp 


Io'y-Mt))dt 


0 


K< 


(2.36) 


The final result, Eq. (2.35) is consistent with the renormalization group with a properly 
defined anomalous dimension matrix T. In particular, the t-dependent 2-vectors of Eq. 
(2.35) are the right-eigenvectors of T, so that we directly obtain, with simple algebra, the 
following conclusions: 


1 . The off-diagonal entry 7 ^^ is given by |I 0 |, p!5 


7» = iVniJ?2(7+)=iV/7l(^)ff2(7l(^)) + NNL 


with H 2 {pi) given by Eq. (2.26). 


(2.37) 
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2. The larger eigenvalue of T, call it 7 ^,?^ is not exactly the leading solution 7 + of Eq. 
(2.29), but is given by 


lY = ^+Mt)) - ^ log 7+\/-X'(7+) 


(2.38) 


where we have included in the anomalous dimension the NL contribution due to the 
running coupling in the coefficient function H{as)- 


3. The lower eigenvalue 7 _ = —§^ 7 qg^ is just the one loop result, so that also 


^ (795 - 7gg) 


(2.39) 


Since the results (2.37) and (2.39) have been discussed elsewhere 130 , we concen¬ 
trate in the following on the new results on 7 _|_ of Eq. (2.29) and (2.38), obtained by 
evaluating the NL eigenvalue Xi- Note that the combination evolves according 

to given in Eq. (2.39) and that the real dynamical contribution in the latter is the 
one due to 7 +, which will be discussed in the following sections. 


On the other hand, the contribution to (2.38) due to the logarithmic derivative of H 
starts from three loops on, and is a peculiarity of the BFKL equation arising from its 
k-fluctuations, which become large when xl approaches the saturating value xl = i/sQ. 


This point was discussed elsewhere |]^ for frozen as and arises from the saddle-point 
features of Eq. (2.11) in the running case. 


To sum up, we have given in this section a definition of the gluon and quark densities 
based on k-factorization only, i.e., we have defined the so-called Qo-scheme. We have also 
derived from the NL BFKL equation (2.1) the anomalous dimension representation (2.35) 
in the regime (2.13), leading to the expression (2.37)-(2.39) of the anomalous dimensions. 


We have thus set up the framework for deriving resummation formulas from the prop¬ 
erties of the kernel eigenvalue on the basis of the definition of 7 + given in Eq. (2.29). 


^Although the t-dependent o^-singularity of the anomalous dimension may be at a slightly different 
position than LUp{t) it remains true that resummation effects are large when 7/, approaches 

1 / 2 . 


13 










3 Real qq Emission Kernel 


The program of computing the NL kernel involves the lowest order computation of the 
2-particle s-channel discontinuity of the four gluon exchange amplitude, of the one loop 
correction to the corresponding 1-gluon intermediate state, and of the purely virtual two 
loop corrections (Fig. 3). Each of these terms can be further divided in sub-terms having 
different colour structure, according to their dependence on Ca, Cp, Nf. In the following 
we will concentrate on the iVj-dependent part, which will be called the qq contribution. 

Up to now all virtual corrections (to one or no-particle discontinuities) have been 
computed by Fadin and collaborators O, 113, while the evaluation of the real emission 


terms has not been completed. In fact, amplitudes P, ^ and squared matrix elements 
for qq 0 and gluonic |]^ contributions are known in the literature, but the phase-space 
integrations have not yet been performed. 

In the following sections we give the details of our recent computation 0 of the qq 
emission kernel, of its eigenvalues, and of the corresponding anomalous dimension. 


3.1 Massive qq production 

In order to regularize s-channel collinear singularities, we will work with massive quarks 
(the relevant formulas in the massless limit and 4-|-2e: dimensions will be given in the next 
section). As an intermediate step, we shall obtain the high energy coefficient function for 
the QQ production via gluon-gluon fusion. 

We thus start by considering the total cross section for the hadroproduction of a 
(heavy) quark-antiquark pair, which at high energies is dominated by the (Regge) gluon 
fusion process of Fig. 4. In the large center of mass energy limit the perturbative expan¬ 
sion of this cross section contains large terms of order [as(M^) log(s/M^)]"’, which need 
to be resummed to all orders. This resummation is performed by using the k-factorized 
expression for the cross section P] 

— - —^(i^kci^k'(j(k, k', M^, ziZ2s)J^i{zi, k)J^2(^2, k'), (3.1) 

A ^2 
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in which the gluon densities Ti satisfy the BFKL equation (2.1), and a is the (lowest-order 
and gauge-invariant) off-shell continuation of the partonic cross section. 


In the language of and s-moments defined in Sec. 2, we introduce the double Mellin 
transformed version of Eq. (3.1), i.e., 

.,,2 d7i<i72<iw ^ 

M a[ - = / - 72 71,72 3.2 

\ s J ioo -H27r)'*7i72 \ QtJ 


where 


rl roc ^21, / n2 \ O' 

= I ( 3 . 3 ) 

d^kd^k' f k^ywk'^y^ ds fM^Y 

= I y[—) ^(k,k,M,.,.2.) 


For large enough the q-integrals are dominated by the BFKL anomalous dimensions 
(2.31), so that the pro cess-dependent high-energy resummation effects are embodied in 
the coefficient function ^^ 02 ( 71 , 72 ), through the cu-dependence of 71 and 72 . Part of this 
section is devoted to the evaluation of this coefficient function. 

Let us start by defining the kinematics of the process under consideration (Fig. 4). 
First of all we introduce a Sudakov parametrization of the exchanged gluons’ momenta 

F ~ zip^ + k^, ~ -Z2pi^ - k'^, (3.5) 


and of the momentum transfer 


= zixip^ - Z 2 X 2 P 2 + A^, (3.6) 

where pi, p 2 denote the (light-like) momenta of the incoming hadrons. Expressing the 
phase-space in terms of these variables, the moments in Eq. (3.4) take the form 

TT 

j - xy - (k - A)2 - m2) X 

X (5(a;i(l-a;2)z/-(-k' +A)2-M2) ( — j A(k,k', A,xi,X2,i^, M'^). 

(3.7) 
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Here A denotes the (off-shell) squared matrix element for the (Regge) gluon fusion channel 
of the QQ production process, and we have dehned 

z/ = ziZ 2 ( 2 pi ■ P2) = Z1Z2S. (3.8) 


In order to understand the colour and singularity structure of the qq kernel, in view of 
the following phase-space integrations, we hnd it convenient to rewrite the squared matrix 
element of Ref. in the following form: 

^ = al{D, + D2 + D3) (3.9) 


where 


= {2Cp - Ca) 


D2 = Ca 


-1 


+ 


(H + C)- 

k2k'2 


- u) 

1 \ . B^ + C^ 

(1 - 2^1 - X2) - O. .n + 


s V — u — t 

2{B - C) 


+ 


D. = C, 


2k2k'2 

(k^(l - X 2 ) + k'^(l - xi) k ■ k') 
2 (k^(l — X 2 ) + k'^(l — xi) -f k ■ k')^ 


k2k'2s 


i/s k2k'2 


and we have introduced the notation 


1 (1 - Xi)(l - X2)z^ Z/(1-Xi-X2) A-(k-hk') 

JD = --1-:-1- 

2 M2 -1 2s s 

X 1 X 2 Z 2 z/(l-xi-X2) A-(k-Fk') 

~2 ~ -u 2i i ’ 

and the Mandelstam variables of the hard sub-process 


(3.10a) 


(3.10b) 

(3.10c) 


(3.11a) 

(3.11b) 


s = {ki + /c2)^, t = u = {ki — k2 — A)^. 


(3.12) 


We have decomposed the matrix element (3.9) in the three terms (3.10) in such a way 
that the integrals over u in Eq. (3.7) are separately convergent. Furthermore, all the 
s-channel collinear singularities are contained in the term in Eq. (3.10c), since the 
various poles of the type 1/s and 1/s^ in (3.10b) cancel each other. This implies that we 
will be allowed to compute the contribution of Di and D 2 to the kernel directly in the 
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massless limit and four space dimensions, and only the term with its simple structure, 
will need a regularization. 

As a hnal remark, note that the hrst term (~ (2Cf — Ca)) is suppressed by a colour 
factor of 1/C^ with respect to the others. We shall see in the following that this sup¬ 
pression factor is relevant in order to assess the magnitude of NL qq contributions to the 
kernel eigenvalue. 

The algebra involved in the evaluation of the phase-space integrals in Eq. (3.7) is 
cumbersome. So we give here only the final result in the a; = 0 limit, and we report the 
details of the phase-space integrations in Appendix A. Following the decomposition (3.9) 
of the matrix element, we split the moment in Eq. (3.4) in three parts and 

as follows 

i/o(7i,72) = i/(')(7i,72) + ^^^'^71,72) + ^^^'^71,72) (3.13) 

where 


= {2Cf - CA)H^\lin2) 

1 - 7i)S(72,1 - 72) 


r(2-71-72) r( 3 -71-72) 


r( 4 - 27 i- 272) r(6-27i- 272) 


(3.14a) 


(3.14b) 


Hi3) ^ 


as „ 5(71,1 - 7i - 72)5(72,1 - 7i - 72)5(2 - 7i - 72, 2 - 7i - 72) 


TT 


2(1 - 7i - 72) 


X 


X 1- 


5 - 271 - 272 


(3.14c) 


and i 7 “"( 7 i, 72 ) is the abelian contribution (2.26). 


The result (3.13)-(3.14) provides, in 7 -space, the leading order coefficient for the pro¬ 
duction of a (heavy) quark-antiquark pair via gluon-gluon fusion, relevant for QQ pro¬ 
duction in hadron-hadron scattering, and completes the results of Ref. where the Cf 
dependent part was computed. 


The above coefficient (in the massless limit) is related to the NL eigenvalue, repre¬ 
senting the action of the kernel on the leading order powerlike eigenfunctions ~ (k^)'’'”^. 
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The massless limit, obtained as the residue of the pole at 71 + 72 = 1 (Cfr. Eqs. (2.25) 
and (2.26)), is singular. In fact, due to the s-channel collinear singularity of D^, the term 
//(3) 

in Eq. (3.14c) has triple and double poles, which must be cancelled by the virtual 
terms. 


3.2 Massless quark limit 


The triple pole singularity of in Eq. (3.14c) is due to the singular behaviour of the 
99 99 cross section of Eq. (3.4) in the region (k —k')^ = ^ k^ ~ k'^. In fact, 

after phase space integration, the partonic cross section approaches the singular limit 


<j^=o(k,k',M^) = 


a. 




tto — / dxi 

^ Jo 


a;i(l - xi) 


qz 


log 1 + 


a;i(l — Xi)o[^ 

M 2 


~ a. 


CXg 

’Gtt 


log 


q 5 

M^~3 


qz 


( 3 . 15 ) 


which is responsible for the triple pole, of the k2-moments in Eq. (3.13). 


As remarked before, this singularity comes from the last term ZI 3 in Eq. (3.10), while 
the s-channel singularities of D 2 cancel each other after phase-space integration. 


It is interesting to notice that such singular behaviour is responsible for a sizeable 
enhancement of the QQ hadroproduction cross section at high energies when the two 
anomalous dimension 71 and 72 approach the saturating value of 1/2. On this basis we 
expect relevant effects due to the resummation of higher orders in the b-quark cross-section 
at Tevatron, which will be the subject of further investigations. 


From the point of view of the qq kernel instead, we have still to include the virtual 
contributions, which are singular also. Having identihed the singular real emission part 
with the term in Eq. (3.14c), we subtract it off and we perform the massless limit 
on the hnite part, leaving to the next section the combination of with the virtual 
radiative corrections, to yield an overall hnite result. 

In analogy with Eq. (2.25) we then dehne the regular qq contribution to the kernel 
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eigenvalue, as follows 


Xrlg = lim (1 - 7 - 7i) (7,7i) + (7,7i)] = 

_CAasf TT y/ 3i2CF-CA) , , 1 + |7(1-7) \ 

Gtt Vsin 7 r 7 j V (3 _ 27)(1 - 27)(1 + 27) J ’ 

(3.16) 


which is obtained directly from Eqs. (3.14a) and (3.14b). 


For 


^ j. \ \ / 

completeness we report the kernel (3.16) in k-space also. By inverting Eq. (2.4), 
as explained in Appendix B, we are able to express the azimuthal averaged real emission 
regular qq kernel as follows. 


reg /, 1^/2 j 


2Cf — Ca NfCn 
Ca 647rk2 


logi(i-l)+2(i + l 


+2 ^C 2 {p) +log-£i(p)^ 

, 1 ’ ^ 

^ 67r k 2 -k '2 


"> L P \P 
22 

7 V 

'k' 2 ' 


+ 


r 

■ p - - 

p. 


+ 


where 

Ci{p) =^arctan(^), 

vP 

1 

^2{p) {Li2{iy/p) - Li2{-iy/p)) , 

p=^, A;< = min(|k|, |k'|), /i;> = max(|k|, |k'|), 

and Li 2 denotes the Euler Dilogarithm. 


(3.17) 


(3.18) 

(3.19) 

(3.20) 
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4 Complete qq Kernel and Anomalous Dimension 


4.1 Combining with virtual terms 


Let us now consider the contribution to the qq kernel of the singular part in Eq. (3.10c), 
whose 7 -moments yield the term in Eq. (3.14c) with the triple pole. The singular 
behaviour for and —>■ 0 of the cross section (3.15) is cancelled in the total NL kernel 
by the corresponding singularities in the virtual terms. 

In order to cancel the = 0 singularity we need the quark loop contribution to the 
one-gluon s-channel discontinuity, which is simply obtained from the NL reggeon-reggeon- 
gluon vertex computed in Ref. |^. By squaring the corresponding amplitude, summing 
over the polarizations of the emitted gluon and averaging over the azimuth, we get 

ctsNf 


1 m 2 


^log^ 
k2_k'2 ^^k'2 


(4.1) 


bvr \ 

in which /i^ is the renormalization scale. By combining with the singular qq emission part 
(3.14) we get the M^-independent expression 


dsNf 

bvr 


log 






qz 


1 . 

k2 _ k'2 k'2 


(4.2) 


which is still singular for q^ = 0 . 

The same computation can be performed in the massless quark limit and 4 -|- 2^ 
dimensions, to yield the ^-dependent kernel (see Appendix C) 


ir2(l-Ke) /q 2 M 1 - - 1 - ^^2 


asNf 

bvr [e: r(l -|- 2e) 


q^ 


eq^ 


1 

k2 - k'2 


( 4 . 3 ) 


In order to cancel the q^ = 0 singularity we need to combine the above results with 
the purely virtual terms, which were computed in dimensional regularization in Ref. |T4 


This can be done by means of Eq. (4.3), and after the cancellation of the singularities at 
£ = 0, the NL qq contribution to the kernel can be expressed directly in 4 dimensions as 
follows (see Appendix C): 


A"o + — ( 1 H log — ) iLo(k, k') -|- 


+ 


a,N 


f 


bvr 


q- 5\ 1 

k2 3 J q2 


R 


k2 _ k'2 k'2 


+ iL,%(k,kO, (4.4) 
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where is the regular part of the real emission kernel given in Eq. (3.17), and we have 
factorized the qq contribution to the running of ag at the scale in agreement with the 
choice of Sec. 2 . 


We remark that such a choice of the scale of o;* is not strictly consequence of our 
NL calculation, since it Exes only the coefficient of the log(l//i^) term, while the upper 
scale of the logarithm can be changed by changing the scale invariant part of the kernel. 
However this choice allows the clear discussion of factorization properties of the BFKL 
equation provided in Sec. 2 (Eqs. (2.35) and (2.36)), and yields the NL coefficient (2.19), 
induced by the running coupling, which differs from unity only from three loops on. 


Finally, we are now able to compute the finite eigenvalue of the scale invariant kernel 
Ki in Eq. (4.4). By using standard identities of type 


1 


108 777 ^ = 


dy 


and 


k2 _ k'2 ° k'2 


d^k ( k_^y 


'o l/k 2 + (1 - y)k 


/2 


= (M2)-"r"5(7,n-7), 


Jq nk'^ \My (M2 + ^k2) 
we obtain the expression 

ttgNf 


QQ 

Xi 


Gtt 


^(Xo(7) + Xo(7)) - ^Xo(7)- 


Ar 2 


TT 


^ \Sm 7 r 7 


3 cos(7r7) 


1 + 17(1-7) 


(3-27)(1 -27)(1 + 27) 


(4.5a) 


(4.5b) 


( 4 . 6 ) 


which is plotted in Fig. 5 as a function of 7 . 


Note that the expression (4.6) is symmetrical for 7 —1 — 7 , apart from the + 0 ( 7 ) term, 
which is due to our choice of k^ as scale of a.^, instead of a more symmetrical combination 
of k and k'. 


4.2 Eigenvalue properties and anomalous dimension 

From the plot of Fig. 5 it is evident that a cancellation between terms with different 
colour structure is at work. In fact, the last term in the eigenvalue (4.6) is essentially 
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the abelian one of Eq. (2.26), except that it is mnltiplied by the small colonr factor 
1/iV^. This means that the non planar diagrams are the ones responsible for the large 
resnmmation effects, giving large contribntions at high valnes of both and k'^. In fact, 
they are snppressed in the large Nc limit when conpled to glnons in i^i, bnt are not when 
conpled to photons in H 2 (Eq. (2.26)). 


From Eq. (4.6) the part of the shift of the pomeron intercept which is proportional to 
Nf can be estimated. For models with a continnnm spectrnm 0 it is simply obtained 
by setting 7 = 1 / 2 , to yield 


Acup — cup 


asNf 

dvr 


21 og 2 - 


AcUf 

wp 


~ —0.05iVjQ!s, 


1 33 TT^ 

8 log 2 


(4.7) 

(4.8) 


which is jnst a few percent for reasonable valnes oi Nf and Once again, the potentially 
large term ~ is snppressed by a colonr factor. 


In order to compnte the NL correction to the anomalons dimensions, we introdnce the 
’renormalized’ conpling 

Noov 


Ota ^ OCg Ota I 1 T 


cnp 


in the dehning eqnation (2.29), so that it takes the form 


a 


1 = — [Xo(7+) + (Xi(7+) - Xi(l/2))] 


UJ 


and then we expand 7 + in the Xo-term aronnd 7 + = 7 '^(q!s), to get 

x\\1l) - Xf (1/2) 


[7+ - = a. 

1 = asXoi.lL)- 


-X'oilL 


(4.9) 


(4.10) 


(4.11) 

(4.12) 


In this way the pertnrbative expansion in powers of ihias) is free of spnrions singnlarities 
for xl = 1/2, which wonld signal precisely a renormalization of the Pomeron intercept. 

The NL correction (4.11) to the leading eigenvalne of the anomalons dimension matrix 
7 + is plotted in Fig. 6 as a fnnction of hiag)- This correction is small compared to the 
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leading anomalous dimension 7 ^, apart from the region of very small 7 ^ corresponding to 
> tts, where it is needed in order to achieve consistency with hxed order renormalization 
group equation. Therefore, higher order effects are negligible. 


A particular comment is needed on how hxed order perturbation theory is recovered 
in this formalism. In fact, since 7 ^ = asjio + O(as^), the small 7 behaviour of Xi{l) 
Eq. (4.11) determines directly the one and two loop anomalous dimensions: 


7 + - 7 l ~ 


CisNf 

bvr 




1 + 


13 


lOA"^ / uj 


a.. 


+ .. 


(4.13) 


and from Eq. (2.31) we obtain 




asNf 

bvr 


23 a, 

1 H-^ + 

6 w 


(4,14) 


which agrees with the two loop computation in the DIS scheme 
More in general, we can say that if 

7+) = — + Aia„ 7^.^^ = Aaa 


Otg 

’ UJ 


(4.15) 


is the NL expression at two loop order, the eigenvalue Xi should have the structure 


/ \ 4I1 44.2 . . 

Xi ( 7 ) — — 5 - H-h (7 —> 1 — 7 ) + running coupling terms 

7 7 


(4.16) 


just on the basis of the 7 = 0 singularities needed to recover (4.14) and of the basic 
k k' symmetry of the BFKL kernel - apart from the asymmetry induced by factorizing 
out as(k^) in Eq. (4.4). 

The singularities at 7 = 0 occurring in Eq. (4.16) are just the collinear ones coming 
from transverse momenta |k| S> A S> |k'|, while the ones at 7 = 1 come from the 


disordered region (|k'| S> A S> |k|), that we have already emphasized as a possible 
source of large higher order effects. The third possible source of collinear singularities - 
the one in the s-channel - turns out to combine with virtual corrections to give rise to 
running coupling effects added in Eq. (4.16). 

The structure of Eqs. (4.15) and (4.16) is interesting for various reasons. Firstly, 
it appears that the NL terms arise from various collinear singular kernels, as already 
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emphasized in Eqs. (2.33)-(2.34) at one-loop level. This is because the two channel 
collinear factorization is not disentangled yet in the BFKL approach (Cfr. Sec. 2.2). 


Secondly, we can follow the suggestion of Ref. |2^, and use the expression (4.16), 
including the (1 — 7 )-terms, as a crude estimate of Xi for fhe gluonic contribution, where 
a complete calculation is not yet available. In that case we should set 


_ asNf 1_ llNgas , _ asNf /5 13 \ 

^ Gtt 127r ’ ^ Gvr \3 6iV^ / 

with the rough estimate 


(4.17) 


lliVoa, f 1 ^ 1 , 

XUM - x,,^u - + (^ 3 :^ j + 


(4.18) 


The last term in Eq. (4.18) is large, and - if extended to large 7 values - would provide 
a negative and sizeable shift of the Pomeron intercept and important resummation effects, 
unlike what happens for the qq contribution. 

Although we do not think the estimate (4.18) is reliable when approaching 7 = 1/2, it 
provides an indication of the existence of unsuppressed planar contributions in the gluonic 
part. The complete calculation is therefore a rather important goal to achieve. 
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5 Discussion 


We have studied in this paper the anomalous dimension representation for the NL BFKL 
equation and the ensuing resummation formulas, by analysing in detail the qq contribu¬ 
tions to the larger anomalous dimension eigenvalue. 

Our first result (Sec. 2) is that the R.G. representation with running coupling and 
resummed anomalous dimensions is valid at NL level in the regime 



( 5 . 1 ) 


that is, if the variable Q;s(t)logi is not too large. If the condition (5.1) is coupled with 
the saddle point estimate of important uj values we end up with the limitation 



(5.2) 


which provides a parabola-like boundary in the log^, log^ plane. 


Conditions of type (5.1) have been already noticed before as singularities of the 


anomalous dimension and those of the type (5.2) have been known for a while to be 


relevant, with different numerical factors, for the occurrence of higher twist unitarization 
effects [27-29]. Here we just emphasize that, within the regions (5.1) and (5.2) there is 
a well dehned resummation of anomalous dimensions, provided here, which is able to 
describe the QCD evolution in agreement with the NL BFKL equation (Fig 7). 

The renormalization group description in the regime (5.1) holds independently of the 
detailed properties of the hard Pomeron, i.e. of the leading singularity in the a;-plane 
which dominates the small-a:, fixed behaviour. The latter occurs in the coefficient 
function, and is dependent on the soft region behaviour of the running coupling, in par¬ 
ticular on its magnitude and shape around = A^. Only if the scale is large enough, 
the Pomeron decouples, and the coefficient function takes a perturbative form, provided 
the rapidity is not so large to allow diffusion from to A^. 

The simplest way to summarize the above results is to write the resummed R.G. 
representation of the DIS structure functions F 2 and Fl. By using formulas (2.35) for the 
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parton densities in the Qo-scheme we obtain 


F^iQ^) 


1 7+^^2(7+) 


Gi + 


.£a 

Cf 




CXs 

Stt 


\7+^^l(7 

where H 2 and Hi are given by Eq. (2.26) and by the relation 


(5.3) 


»i(7) = 

1 +17(1-7) 

respectively, and G±^ are the R.G. expressions of Eq. (2.36). In particular, 
1 


(5.4) 


G^ = 


: exp 


7+(a^(t))dt 


(5.5) 


7+V-X'(7+) 

contains the perturbative coefficient of Eq. (2.19), which provides an additional contri¬ 
bution to the effective anomalous dimension of Eq. (2.38). 


Note that, once a complete NL computation of 7 + will be available, all the relevant 
coefficient and anomalous dimensions in Eq. (5.3) will be known in the Qo-scheme, so 
as to provide a factorization scheme independent expression for the measurable structure 
functions. 


In fact, relating Eqs. (5.3) and (2.35) is equivalent to computing the quark and 
gluon coefficient functions in the Qo"Scheme to all leading orders. By simple algebra, and 
neglecting subleading contributions, we get 


Gt 



1 , 


GpNf 

Ga 


't+HAiA 



7 = o; (5-6) 

7 = 7+^^L(7+). (5.7) 


In this paper we have further provided a computation of the gg-contribution to the 
BFKL kernel, proving the running coupling factorization, and we have evaluated the 
corresponding NL resummation in the anomalous dimension eigenvalue 7 +. 

We hnd that higher order effects for the Nf dependent part are small in 7 _|_, while 
they are not small in the coefficients H 2 and H^. This difference is due to the non-planar 
nature of the diagrams which yield the most important large k contributions. They are 
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suppressed by a colour factor when coupling to gluons, while they are not when coupling 
to photons. 


We have also found that running coupling effects, although expected, are particularly 
important. First of all, they affect the coefficient functions by factors of type H in Eqs. 
(2.19) and (5.5) due to the fluctuations of the anomalous dimension variable 7 , which are 
large when 7 approaches the saturation value 7 = 1/2. Ultimately, this is the reason for 
the breaking of the R. G. representation itself when approaching the critical value (5.1), 
as also noticed in a recent paper . 


Furthermore, the running of dgit) ~ 1/ht at large t values, emphasizes the diffusion 
towards small values of k^, and the need of smoothing out the effective coupling around 
the Landau pole. This in turn clarihes the fact that the bare hard Pomeron singularity is 
actually dependent on soft physics, even if small-x scaling violations are not. 


It thus appears that, so far, apart from running coupling effects, large higher order 
contributions only occur in the jqg entry of the anomalous dimension matrix, which was 
the basis for an early explanation of the HERA data 0]. Note, however, that the absence 
of higher order effects in 7 _|_ may be a feature of the A^/-dependent part only, which is 
mostly nonplanar. For the gluonic part, planar diagrams could also contribute, as it 
happens in the crude collinear estimate of Eq. (4.18). Such possibly large effects renew 
the present interest in a complete computation, which is hopefully to be obtained soon. 
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Appendix A: Kinematics and Integrals 


We give here some details of the computation of the Hq function (3.13). 

Starting from the Sudakov parametrization (3.5)-(3.6) of the relevant momenta, we 
write the two-body phase space of the produced qq pair as follows 


2 ) = — -dxidx2d'^AS{ x 2{1 — xi)u — (k — A)^ — M^) x 

X (5(a;i(l — X2)i' — (—k' -|- A)^ — M^) = 

1 dxi d^A^(_ aVm2 _ ^\ _ 

Stt a;i(l — xi) 27r y xi{l — xi) ^ j 

1 dx2 ^ + M'^ 2 ^ 

Stt 0:2(1 - 0:2) 27 r y 0:2(1-0:2) ^ j’ 


(A. la) 
(A.lb) 

(A.lc) 


where 


A = A — kxi — k'(l — 0 : 1 ), (A.2a) 

A = A — k'o :2 — k(l — 0 : 2 ).. (A.2b) 


So, since the hard cross section is related to the matrix element by 

(j(k, k',M^,z/) = ^ J d^{l, 2 )A(k,k',A,xi,X 2 ,M‘^), (A.3) 

we get, by (A.la), the expression (3.7). 

Since both parametrizations of the phase-space (A.lb) and (A.lc) are needed in the 
computation of the integrals, for each of them we give the expression of the relevant scalar 
quantities. Starting with the former, we have 


^ A 

0 : 1(1 - 0 : 1 )’ 

(A.4) 

m2 - i =(1 - 0 : 1 )(s + k2) + o:ik'2 - 2A ■ k'. 

(A.5) 

— u = 0 : 1(5 -I- k)2 -I- (1 — o:i)k'2 — 2A ■ k'. 

(A.6) 

r ~ 12 

(1 —o:i)q —A -I- m2 


(1 - Xi)v 

(A.7) 
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for the latter 


. _ A% M2 
^ 0:2(1-xa)’ 

M^ — i =(1 — a:2)(s + k'^) + — 2A ■ k, 

M^ — u =X 2 {s + k'2) + (1 — a: 2 )k 2 + 2A ■ k, 


Xi 


(1 - a: 2 )q + A 


+ M 2 


(1 - X2)iy 


and for both of them 


z/ =s + q 2 , 

q =k - k'. 


(A.8) 

(A.9) 

(A.IO) 

(A.ll) 


(A.12) 

(A.13) 


In order to compute Hq, we find it convenient to choose, for each term in the matrix 
element (3.9), the most suitable phase space parametrization of type (A.lb) or (A.lc). 
The choice is done by requiring that the denominator of each term in the matrix element 
(3.9) contains at most two of the invariants given before, and that at most one non-trivial 
azimuthal integration is needed. It is straightforward to check that at least one of our 
parametrizations satishes these requirements for every term in Eqs. (3.9) and (3.10). 


With these premises we hrst use the ^-function in (A.lb) or (A.lc) to eliminate the 
variable z/, then, with the help of a Feynman parametrization of the denominators, we 
perform the integration over the transverse components of the momentum transfer A. In 
order to compute the Hq function, the best way is to evaluate the k^-moments at this 
stage, since all the integrals needed can be reduced to the forms 

r^d^k/k^y 1 


0 nk^ [My (M^ + ^ky 




(A. 14) 


k 2 


72 


/g Trkf 7rk| \M‘^ 

_ ,,^ 2 ^-n^- 7 l^- 72 ^( 7 l)^( 72 )^(r^- 7 l -72) 


My {M^ + ^,kl + ^2kiy 




r(n) 


(A.15) 
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71 


k 2 

M2 


72 


(A.16) 


d^ki d^k2 / k? \ 
Jq Trkf nkl \M‘^) 


1 

[M2 + e(ki + k2)2]- 


. 2 N-n^- 7 i- 72 ^(^ + ^ ~^2)r(i -7i -72)r(7i)r(72) 

^ ^ ^ r(n)r(i-7i)r(i-72) 


The remaining integrals over the longitudinal momentum fractions and the Feynman 
parameters are evaluated in a straightforward way in terms of Gamma and Beta functions, 
the only trouble being the number of terms, about one hundred. After doing the Gamma 
-function algebra the result (3.13) is found. 
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Appendix B: The kernel in k-space 


The inverse Mellin transform needed to compute the regular kernel (3.15) in transverse 


momentum space involves integrals of the type 



- X ' 

2tti 



(B.l) 


where x = k^/k'^ and R{'y) is a holomorphic function of 7 , symmetric under the trans¬ 
formation 7 —> 1 — 7 , which is provided by the eigenvalue (3.16). 


Integrals of this type can be evaluated by closing the contour on the 1. h. poles for 
a; < 1 and on the r. h. ones for a; > 1, so that they reduce to inhnite series. Due to the 
simple structure of all the integrals can be expressed in terms of the following basic 
series 


00 


n=0 


1 

1 — a;’ 


00 


E 

72=0 


(_l)na,n 

n+ i 


00 


E 

72=0 


(_l)na,n 



arctan(^/a;), 


2 

i^/x 


[Li2{iy/x) - Li2{-iy/x)) , 


where Li 2 denotes the Euler Dilogarithm. 


(B.2) 

(B.3) 

(B.4) 


By introducing the expression of the regular eigenvalue (3.16) in (B.l) and performing 
the 7 -integrals as described above, the expression (3.17) for the kernel in k space follows. 
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Appendix C: Dimensional Regularization 


The peculiarity of working in 4 + 2e dimensions is that the kernel is dependent on e as a 
regularization parameter, while its eigenvalues are hnite (for e —0). Therefore, in order 
to compute the qq contribution to the kernel eigenvalue in dimensional regularization we 
have to compute the action of the (e-dependent) real and virtual contributions on a set 
of test functions. By adding the various terms, the cancellation of collinear s-channel 
singularities comes out, and the hnite qq contribution to the NL eigenvalue is obtained. 
Since we want to show that the NL kernel has the form (2.1), we hnd it convenient to 
apply the kernel to the powerlike eigenfunctions of the leading kernel, i. e. 

/.(k^) = (k^)^-^ (C.l) 


C.l Real qq emission 

Let us hrst notice that the collinear singularities are all contained in the term in 
Eq. (3.10c), so that Di and D 2 can be computed directly in 4 space dimensions by taking 
the residue at the simple pole for 71 -|- 72 = 1, along the lines of Eq. (2.25). The ZI 3 term 
can be evaluated as follows. 

After a simple algebra we hnd 

Ds = 4 [1 - 2(1 - xif + 2 x 1(1 - xi)] + (C.2) 

1 / 0 '- 


+T or — terms, (C.3) 

s S"* 


where in the last line we refer to terms of the form 
P(Xi)(k2)“(k'2)/3 

gn 


(C.4) 


(where P(xi) is a polynomial in xi), which just vanish in dimensional regularization. 


By multiplying Eq. (C.2) by the eigenfunction of the real kernel in Eq. (C.l) and 
introducing the 4 -|- 2e dimensional phase-space integrals over the internal momentum 
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transfer A and the transverse momentnm of the emitted qq pair, we get 


Ds(^f^{k) = 

= (47r/i^)“^^asiV/ J j dxi{l — Xi)dy[(k + x 

l + 2xi(l-a;i)-2(l-a;i)^ 

[A2 + a;i(l-a;i)|/q2p ’ 

where a Feynman parametrization of the denominator ~ l/{us) has also been introdnced. 
The integrals in (C.4) can be evalnated by hrst integrating over the two transverse mo¬ 
menta A and q, and then on the longitudinal momentum transfer Xi and the Feynman 
parameter y. The hnal result is easily expressed in terms of F-functions, as follows 


D,®f^{k) = 


/ k 2 Y'V/a, 1 F(1-7-2£)F(7-F£)F(1 + 2 £)F 2 ( 1 +£) 

1 

2 (1+£) 1 

\ATiy?) IbTT F(1 — 7 )F( 7 -|-3^) 

_F(2 + 2£) 

F(4 + 2£)_ 


(C. 6 ) 


and has the following expansion around £ = 0 

2 e 


Airy? 
10 


Nfas 


127r 


+ ^ (^2x0(7) “ - 7) - 4 V''( 7 ) + 


+ 2 X 0 ( 7 ) - ^Xo(7) + 


28 


(C.7) 


Note that in 4 -|- 2^ dimensions, due to the lack of scale-invariance of the kernel, the 
powers of are no longer eigenfunctions of the kernel, and are in fact multiplied by 
powers of Anyway, after the combination with virtual terms, the £ —0 limit 

will be taken, and it will be possible to isolate the non scale-invariant terms (related to 
the running of ag), from the scale-invariant NL corrections. In performing this procedure 
we are obliged to expand all partial results up to relative order e^, because of the double 
poles occurring in Eq. (C.7). 

C.2 Vertex corrections 

The vertex corrections to the 1-gluon production amplitude have been computed in 
Ref. [j^ (Vj-dependent part) and |]^ (Vj-independent part). In the following only the 
^/-dependent part will be considered. 
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By squaring the 1-gluon production amplitude, together with its one-loop corrections 
and summing over the emitted gluon’s polarization, the contribution to the kernel in 


4 -h 2£ dimensions is readily obtained, as follows^ 


NfOs 

Git yAnfi 


{4n)-^T\l-e) 


1 1 
£ q 2 


1 

k 2 - k '2 k '2 


+ 


terms of order £ and hnite for q = 0. 


(C.8) 


The action of this part of the (regularized) kernel on the leading order eigenfunctions 
f-y(k) is obtained by performing the integration over q in 4 -|- 2e dimensions. However, 
only the hrst term in (C.8) needs be regularized, while the second one is hnite in both 
the £ = 0 and q = 0 limit. 

The second term is easily computed in four dimensions, by using the integral repre¬ 
sentation 


1 

k2 - k'2 k'2 


dy 


(C.9) 


/o + (1 - 2/)k'2 

and then performing the k' and y integrations according to the standard techniques of 
Appendix A, to yield the result 

^(^)'. (C.IO) 

DVT \Sm 7 r 7 / 

The hrst term in Eq. (C.6) is computed by means of a Feynman parametrization to 
combine the two powers 1/q^ and [(k-f q)^]''"^, and then performing the 2-|-2£-dimensional 
integration over q with the standard techniques of dimensional regularization. The hnal 
integral over the Feynman parameter is straightforward so that the contribution of the 
kernel (C.6) can be expressed as 


1 r(l - 7 - £)r(l + £)r(7 -h e) agNf 


£2 r(i-7)r(7 + 2£)r(i-£) 

'asNf \ ( TT ^ ^ 


Gtt \47r/i' 


r2(l-£)(47r)-^ + 


Gtt 


sm TT'y 


(C.ll) 


^Following Ref. ||^, we define the renormalized charge in the MS renormalization scheme as follows, 
9 9f^y V 3 (47r)2+2e2e 
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In the £ = 0 limit, (C.ll) has the expansion 


asNf 

hvr 





4 + ^ (xo(7) + ^^'( 7 ) 


£■ 


£ 



(C.12) 


C.3 Trajectory correction and finiteness of the kernel: 


The purely virtual two loop contribution was computed in Ref. 
dependent part has the following series expansion around e = 0 


M 


<52+2^ (q) 




Nfas 

Utt 


'2(47r)- 


+ 


dvr/i^ 


1 51 

—^ + 


28 


3£ 9 


The Nf- 

(C.13) 


By integrating the above result over the phase-space and by combining it with 

the real emission contribution (C.7) and the vertex correction (C.12), the action of the 
total NL qq kernel on the leading order eigenfunction f-y{k.) is finally found, as follows 




NfOs 

Gtt 



5 

3 


Xo( 7 ) + 


— f -^] b cos(7r7) ^ + 


-f 0 (e). 


(C.14) 


Note that only the first term is not scale-invariant, and it is proportional to the lowest 
order eigenvalue, with a coefficient which is just the qq contribution to the R. G. beta- 
function. 

Finally, by combining the real emission terms in Eq. (C.7) with the virtual ones in 
Eqs. (C.IO), (C.12) and (C.14), the results in Eqs. (4.4) and (4.6) of the text are obtained. 
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Figure 1: BFKL equation for the gluon density up to next-to-leading level. Wavy lines 
denote (Regge) gluon exchanges. 



Figure 2: (a) Low energy and (b) high-energy part of the quark density in the DIS-Qo 
scheme. Dotted lines denote the electroweak probe and an F 2 projection is understood. 
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Figure 3: (a) qq and (b) low energy gg contributions to the NL kernel, together with (c) 
one-loop corrections to Ig state and (d) two-loop virtual corrections. 
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Figure 4: qq contribution to the 4 (Regge) gluon amplitude absorptive part. 
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Log(QV) 


Figure 7: Qualitative plot of the anomalous dimension regimes in the \ogx/\ogQ^ plane. 
Curve I is the boundary for the validity of resummation formulas. 
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